Thermodynamically Stable Blue Phases 
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We show theoretically that flexoelectricity stabilizes blue phases in chiral hquid crystals. Induced 
internal polarization reduces the elastic energy cost of splay and bend deformations surrounding 
singular lines in the director field. The energy of regions of double twist is unchanged. This in turn 
reduces the free energy of the blue phase with respect to that of the chiral nematic phase, leading 
to stability over a wider temperature range. The theory explains the discovery of large temperature 
range blue phases in highly flexoelectric "bimesogenic" and "bent-core" materials, and predicts how 
this range may be increased further. 
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The blue phases (BPs) arc liquid crystalline; they have 
orientational molecular order, yet flow like liquids. Up to 
three thermodynamically distinct BPs are observed upon 
cooling from the isotropic liquid: BPIII, BPII, and BPI 
respectively. BPIII is amorphous, in that it lacks long- 
range translational order, whereas BPII and BPI are pe- 
riodic in three dimensions. They exhibit a frustrated 
structure consisting of a network of singular "disclina- 
tion" lines in the field describing the average orientation 
of the molecules - the "director" field. The periodicity is 
typically of the order of the wavelength of visible light. 
This leads to vivid colored Bragg-like reflections, and a 
partial 3D photonic band gap. They are thus perhaps 
the only example of self-assembled 3D photonic crystals 
which, because of their fluidity, have optical properties 
that are readily switchable in an applied electric fleld 

Potential device applications include 3D lasers [1,T3| 
and displays 0. 

Temperature stability is crucial to many applications; 
yet, the width of the entire BP region is only ^ 1°C in 
typical liquid crystal (LC) materials. Polymer stabiliza- 
tion has been employed to widen this temperature range 
0; however, the polymer network leads to restricted tun- 
ability and does not represent a true thermodynamic sta- 
bilization of the phase. BPs stable over a range of up to 
50°C, including room temperature, were reported in 2005 
[loj . The realization of polymer stabilized and nonpoly- 
mer stabilized BPs with wide temperature ranges has led 
to a resurgence of interest in the area. In the latter case a 
new theoretical challenge arose: to understand the mech- 
anism of stabilization. In this Letter we demonstrate that 
it is a result of internal flexoelectric polarization. 

Flexoelectricity in LCs is a linear coupling between 
applied electric fleld and induced distortion [llj . The in- 
verse effect results in naturally distorted director struc- 
tures, such as the BPs or LC colloids, having a flex- 
oelecrically induced internal polarization field. This 
deformation-induced polarization has the effect of renor- 



malizing the curvature elastic energy of director distor- 
tions. Recently, materials with nonsymmctric molecules 
have been developed with an unusually lar ge fl exoelectric 
response [l2- 14|. The stable BPs of Ref. [1^ were com- 
posed of "bimesogenic" LCs, for which very large flexo- 
electric coefficients have been reported 12 1. 



Upon cooling, the BPs usually undergo a transition 
to the helical cholcsteric, or "chiral nematic" phase (de- 
noted N*). The stability of the BPs is primarily a mea- 
sure of their free energy with respect to that of the N*. 
The fact that the BPs can be stable compared to the 
defect- free N*, even in a narrow temperature range, is 
initially somewhat surprising: they contain a dense net- 
work of disclination lines that are energetically unfavor- 
able due to hig h elastic distortion. Following Meiboom 
et al. 15|, llGj the stability can be understood by the 



fact that the N* is not locally the lowest free energy 
conflguration for chiral molecules. In the N*, the chi- 
ral molecules reduce their free energy by twisting with 
respect to each other along the helical axis. Fig. [Ija). 
However, molecules within pseudo planes perpendicular 
to that axis are still constrained to lie such that their av- 
erage orientation is parallel to each other. A lower energy 
conflguration is achievable in regions of "double twist". 
Such regions cannot extend indeflnitely without topolog- 
ical defects, though they can extend parallel to the direc- 
tor, and a short way in directions perpendicular to this, 
to form cylinders of double twist, Figs. [IJb,c). BPII and 
BPI may be considered to be composed of these struc- 
tures, called double twist cylinders (DTCs), arranged in 
regular cubic lattices, e.g., Fig.md). The DTCs have free 
energy lower than the corresponding N*. In between the 
DTCs are a network of disclination lines of topological 
charge s = —1/2. Because of the energy cost of elastic 
distortions, the free energy in the vicinity is higher than 
that of the corresponding N*. If the energy saving from 
the DTCs outweighs the energy cost due to the disclina- 
tion lines at a given temperature, the BP will be stable 



2 




FIG. 1: (Color online) (a) Single twist in the chiral nematic 
phase, (b) Regions of double twist may form double twist 
cylinders (DTCs). (c) Cross section view of a DTC. (d) BPI 
is composed of a cubic arrangement of DTCs with space group 
0*"(/4i32). 
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(a) Electric/Polarization field 



FIG. 2: (Color online) (a) Flexoelectric materials may form a 
characteristic splay-bend deformation of the director n(x) in 
an external electric field [HI. [l7|. (b) The s = —1/2 disclina- 
tion lines are regions of alternating splay and bend. 



with respect to the N* at that temperature [15|, ll6| . 

It is interesting to note that the director profile of the 
s = —1/2 disclination is the splay-bcnd deformation, 
Fig. [2{b), characteristic of flexoelectric materials in an 
external electric field [HI, [l3|: rotated onto itself, Fig. [5] 
Thus, one may intuitively expect that highly flexoelec- 
tric materials would form such a distorted structure more 
readily. Yet, in the absence of external electric fields, a 
theory based on the local free energy density must be 
completely determined by the elastic coefficients. In this 
case, the stability of the BPs must be the result of a par- 
ticularly favorable combination of clastic coefficients. By 
explicitly considering the energy of the internal polariza- 
tion field, its effect on renormalizing (in effect, reducing) 
the values of the corresponding elastic coefficients may be 
investigated [THIil. To second order in spatial gradients 
of n, and second order in induced electric displacement 
D, the free energy density may be written in symbolic 
matrix notation: 



/ = k?t?n+lk°(an) 



(1) 



where tensors and k^* arc clastic coefficients at con- 
stant D, corresponding to the linear (chiral) and square 
of director gradients, e is the local dielectric permit- 
tivity tensor, and e represents the flexoelectric coeffi- 



cients, which couple polarization to distortion. The di- 
mensionality of boldface symbols is to be assumed, e.g, 
k^ (c?n)^ may be written in a specific coordinate frame as 
diUjOkni, where derivatives are with respect to po- 



^2 ijkl 



sition, and we use summation over repeated indices. The 
equilibrium value of D is found by minimization which, 
upon substitution back into /, gives 

/ = k?c?n+l(k°-e2e-i)(c?n)^ (2) 

Imposing the locally uniaxial chiral symmetry of the N*, 
characterized by the point group Doo , and using the fact 
that n is a unit vector, / may be written in the more 
traditional form [13, HH 



1 2 

-K2 (n • V X n + gr) 



(n X V X n) 



+ - {K2 + i^4) V • [(n • V) n - n (V ■ n)] . 

(3) 

Ki, K2, and K3 are the splay, twist, and bend Frank clas- 
tic coefficients, the reciprocal length scale q determines 
the local chiral twisting power, and e|| and are the 
corresponding components of the dielectric permittivity. 
The combination (A'2 + K4) is the "saddle-splay" elastic 
coefficient ei and 63 are the flexoelectric coefficients 
as defined by Meyer [ll|. We see that, for nonzero flex- 
oelectric coefficients, Ki and K3 are reduced, while the 
twist coefficient K2, and K4 are unchanged (which is the 
consequence of splay and bend flexoelectric symmetry). 

The reduction of Frank coefficients is ~ 1% [3l in 
typical LCs, and therefore usually ignored. However, it 
will be more significant for the highly flexoelectric LCs 
now available (e.g. ~ 20%, see below). The effect of this 
on the stability of the BPs may be understood as follows: 
the "renormalization" of Ki and K3 leads to a reduction 
of the energy cost of the disclination lines. Since K2 
and K4 are unchanged, the energy saving due to double 
twist is unchanged. It has previously been shown, using 
numerical calculations, that the range of stability of the 
BPs is reduced as the ratio {Ki + K^)/K2 is increased 
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In order to derive an analytic model to quantify the 
flexoelectric BP stabilization, we assume the following 
simplifications: the magnitude of the flexoelectric coeffi- 
cients are equal |ei| — |e2| = e, the LC is dielcctrically 
isotropic ey = e±^ = e, the unrenormalized elastic coef- 
ficients are equal iff = K2 = KJ^ = K, and K4 = 0. 
None of these assumptions alter the qualitative effect, 
but of course for a full quantitative agreement one would 
need to take into account particular different values of 
these coefficients in particular materials. Further, we ig- 
nore the effects of induced space charge, and nonlocal 



3 



electric field effects [2^. The effective splay and bend 
elastic coefficients are now equal to rK where r is the 
single renormalization factor from Eq. ([3]) 



1 - eV(A'e) 



(4) 



Following Meibooni et al. [l5|, |25[ , the stability of the 
BPs may be investigated by considering the energy of a 
single s = —1/2 disclination line, required to topologi- 
cally match neighboring regions of DTCs. The strength 
of this approach is that it leads to an analytic expres- 
sion for the temperature range of stability, though there 
may be limitations due to the inherent simplifying as- 
sumptions. The director field is approximately given by 
n = (cos (s(/)) , sin (s0) , 0), where 4> is the polar angle 
in local cylindrical coordinates. Three contributions to 
the free energy per unit length are required. (1) The 
thermodynamic free energy cost of the core of the de- 
fect. This is proportional to the cross-sectional area of 
the core i^core = airR^, where R is the core radius, a 
is the energy density of the core, which is a function of 
the actual temperature relative to the isotropic transition 
temperature. (2) The elastic energy cost of splay and 
bend deformations surrounding the disclination, which is 
dependent on the renormalized coefficients Ki and A'3. 
Including the flcxoelectric renormalization, this energy is 
Fsh = (l/4)7rr A' In {Rnmx/R) per unit length of disclina- 
tion. Here, i?max is the outer radius of the disclination 
line and is of the order of the distance between disclina- 
tion lines in the BP. (3) The energy saving due to double 
twist. This is encoded through the volume integral of the 
splay-bend term in Eq. (jS]), converted to a surface inte- 
gral over the surface of the core. Note that this term is 
dependent on the unrcnormalized Frank coefficients K2 
and Ar4, and can be shown in the one-constant approxi- 
mation to be given by -Ftwist = pcr unit length of 



a DTC [151 . The total energy F^, 



■ twist 



+ is 



1 



1 



F = airR^ - -ttK + -BnrK In (i?,„ax/i?) • (5) 

A parameter B is included to compensate for the inaccu- 
racy of the analytic expression. Numerical calculations 
lit 25 1 have estimated that B « 0.5, a value that we 
now adopt. As a first approximation, the coefRcient a 
is linear in temperature T, i.e., a — a{Ti — T) = a AT 
for T < Ti , where Ti is the isotropic transition tempera- 
ture and a may be estimated from the latent heat of the 
isotropic-N* transition [l^. The BP will be stable with 
respect to the N* when F is negative. The core radius 
R is determined by minimization of F, giving the ratio 
Rc ~ \\JrK / a. which is essentially the nematic corre- 
lation length, reduced by the renormalization factor r. 
Substituting this back into the free energy, the range of 
stability becomes 



where A = K/ (l6ai?^^^) . 

It is instructive to calculate the temperature range pre- 
dicted by Eq. ^ for typical, experimentally determined, 
material parameters. Care must be taken: do the exper- 
imental values represent the "bare" or the renormalized 
elastic coefficients? In general, the values of the elas- 
tic coefficients will depend on the conditions of measure- 
ment, and be related [261 according to 



= X--ef/e(,=AT-ef/x|| 



(7) 



and Ki 



Superscripts E, D, and P denote measurements at con- 
stant electric field, displacement, and polarization respec- 
tively. Most experiments that measure these parameters 
employ the application of an alternating electric field, 
in which case constant- £) coefficients are measured [s^ . 
The ratio |ei — eaj/AT can be measured via "flexoelectro- 
optic effect" experiments, also under the application of 
high-frequency electric fields [13, 121- Since ei and 63 



must generally be of opposite sign for a periodic splay- 
bend pattern to be formed [ll| . we assume this results 



in the ratio 2e/K in our approximation 33|. When com- 
bined with Freedericksz transition [2^ measurements, e, 
K and e may be individually determined. 

Consider the typical values a = 8 kJ m"-^, i?max = 
100 nm 

and experimental values for the bimeso genic 
material "FFOllOCB" (|e| = 10.6 pCm-\ K = {Ki + 
K^)/2 = 8 pN, e = 7.75 £0, measured in [Hi). These give 
a 20% flcxoelectric reduction of ATi and (r = 0.8). 
Equation ^ then predicts a temperature range of sta- 
bility of AT = 43 K, which is an order of magnitude 
greater than in usual chiral LCs, and is in approximate 
agreement with experiment on mixtures of such materi- 
als [l^l- For comparison, if flcxoelcctricity is not taken 
into account in the above analysis, the predicted range is 
AT = 7 K. Now consider the standard LC "7CB," with 
experimental values measured using the same methods 
del = 2.4pCm-\ A:= (fs:i-hA:3)/2 = 4.9 pN, e = 10.5eo 
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AT = ^rcxp[(8/r) - 1] 



(6) 



28|). This gives a flcxoelectric renormalization factor 
r=0.99, and the predicted range of BP stability increases 
from AT = 4.2 K to only AT = 4.6 K. Again, this is 
comparable with experimental results. 

The model also explains two observations in bent-core 
LCs: (1) Blue phases may be induced by doping with 
bent-core molecules 2^. (2) Unusually stable BPs exist 
in bent-core LCs [3C|. This may now be understood by 
the fact that, like bimesogens, bent-core LCs have un- 
usually large flcxoelectric properties 13 1. 

As r is reduced further, AT can become very large 
(for example, larger than the transition temperature Ti 
itself). This is not unphysical, rather it indicates that 
the BP is stable with respect to the N* at all tempera- 
tures below the transition from the isotropic. In this case 
the BP will transition directly to the lower-temperature 
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smectic or crystalline phase, at which point the theory 
is cut off. The model may be refined to account for the 
nonlinear dependence of the thermodynamic free energy 
density of the disclination core on temperature far from 
the transition. The core energy a measures the differ- 
ence between the ordered bulk and the disordered core. 
It scales with temperature as the (square of) the or- 
der parameter does, and at low temperatures becomes 
constant. To take this into account, we may instead 
let a = aTi (1 — exp [T/Ti — 1]) (or any another func- 
tion that saturates in the desired manner). In this case 
the predicted range of stability is greater still (46 K for 
FFOllOCB using the above values). 

A numerical minimization of the free energy, such as 
that undertaken in Refs. but which further in- 

cludes the flexoelectric renormalization, would be illu- 
minating. In this way, the assumptions used to gener- 
ate the analytic theory could be relaxed. In particular, 
it should be noted that our analysis of the disclination 
core is highly simplified, and the accuracy of the predic- 
tions may be improved by using a tensor order param- 
eter that will account for the specific symmetry of the 
core. Nevertheless, our analytic model makes clear the 
strong dependence of the temperature range of stabil- 
ity on the flexoelectric coefficients, and the inherent ap- 
proximations are somewhat validated through the above 
comparison with experiment. The theory makes two key 
predictions: (1) Further stabilization of the BPs can 
be achieved by "molecular engineering" of more highly 
flexoelectric LCs. For example, by further exploiting 
the known relationships between molecular structure and 
flexoclectricity (sec, e.g., Refs. dMUll). (2) The depen- 
dence of the splay and bend elastic coefficients on the 
flexoelectric coefficients is large in certain materials, and 
should be accessible to experimental verification. Indeed, 
for highly flexoelectric materials, the flexoelectric proper- 
ties may be accurately determined simply by measuring 
the elastic and dielectric properties under different con- 
ditions, according to Eqs. ((T)). 

We conclude that ficxoelectricity strongly affects the 
energetics of distorted equilibrium structures. The in- 
ternal polarization field affects all splayed and/or bent 
structures (since all LCs are, to some extent, ficxoelec- 
tric), but is particularly relevant in highly flexoelectric 
LCs made of nonsymmetric molecules. This accounts for 
the existence of BPs that are stable over an unusually 
large temperature range. 
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